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We consider non-equilibrium time evolution after a quench of a global Hamiltonian parameter in 
systems described by Hamiltonians with local interactions. Within this background, we propose a 
protocol that allows to change global properties of the state by flipping a switch that modifies a local 
term of the Hamiltonian (creating a defect). A light-cone that separates two globally different regions 
originates from the switch. The expectation values of macroscopic observables, that is to say local 
observables that are spatially averaged within a subsystem, slowly approach new asymptotic values 
determined by the defect. The process is almost reversible: flipping again the switch produces a 
new light-cone with the two regions inverted. Finally, we test the protocol under repeated projective 
measurements. As explicit example we study the dynamics in a simple exactly solvable model but 
analogues descriptions apply also to generic models. 


In quantum mechanics the issue of locality/causality is 
very delicate and gave rise to brilliant debates extending 
over a century of quantum physics^. However in some 
contexts causality does emerge clearly. A shining exam¬ 
ple is given by the Lieb-Robinson bound^, which states 
the existence of a maximal velocity at which information 
propagates in many-body systems with sufficiently fast 
decaying interactions. The typical manifestation is the 
emergence of light-cones in the time evolution of local 
observables^. Importantly, over the time some hypothe¬ 
ses behind the original proof have been made milder^ and 
light-cones have been observed also in situations where, a 
priori, a similar behavior would have not been expected®. 

This kind of locality in quantum mechanics is solidly at 
the base of the descriptions of subsystems at late time af¬ 
ter global quenches in many-body quantum systems. The 
current understanding is that nonlocal conserved quanti¬ 
ties do not affect the late time (stationary) expectation 
values (EVs) of local observables after a global Hamilto¬ 
nian parameter is changed: the effective ensembles de¬ 
scribing the late time physics can be obtained by maxi¬ 
mizing the entropy under the constraints of the local and 
quasilocal conservation laws® (in quantum field theories 
the construction is more delicate'). This picture has the 
remarkable consequence that globally different ensembles 
describe the stationary properties of local observables 
evolving with Hamiltonians that look very similar to each 
other. A small integrability breaking perturbation is in¬ 
deed sufficient to reduce the infinite set of (quasi)local 
charges of the integrable model to a finite set, which in 
the most generic situation consists of just the Hamilto¬ 
nian. The EVs of local observables experience in turn 
a crossover^®^^® from a quasi-stationary behavior, which 
can be captured by a generalized Gibbs ensemble®, to the 
actual stationary value, describable instead by a Gibbs 
ensemble®. These are usually known as prethermaliza- 
tion^®^^® or pre-relaxation®®’®! behaviors. We stress that 
this picture applies also to the case where the unper¬ 
turbed model is generic but there is at least a local charge 
besides the Hamiltonian. 

It is well known that also local perturbations can break 


global symmetries of a generic model or even integrabil¬ 
ity (in fact, a standard problem is to identify the local 
perturbations that preserve integrability®®). 

A (quasi)local conservation law of the unperturbed 
model (which can be either integrable or generic) can 
react to a local defect in two principal ways: (i) it main¬ 
tains its bulk part unchanged, (ii) it becomes extinct. 

The first aim of this paper is to point out that the 
extinction of local charges results in a crossover between 
globally different states. Assuming local relaxation, this 
can be deduced from the time evolution of the EV of a 
charge Q that becomes extinct. It turns out that the time 
derivative of the EV of Q can not be recast as a commu¬ 
tator between the Hamiltonian and some localized oper¬ 
ator. Thus, generally the EV of the commutator remains 
nonzero also at late times and, in turn, the EV of Q 
displays a typical linear behavior. One can then exhibit 
observables (essentially, the current associated with Q) 
with different EVs inside and outside a light-cone, even 
far away from the defect (see Sec. I): the light-cone sep¬ 
arates regions with globally different properties. 

We use this feature to engineer a local switch that al¬ 
lows one to change the state globally in an almost re¬ 
versible way. The protocol is the following: 

1. The system is prepared in the ground state I'E'o) of 
a translation invariant Hamiltonian Hq. 

2. At time tq a global parameter is changed and the 
state evolves with a new Hamiltonian iJoff (global 
quench). Let us call V a local perturbation that 
spoils a set of local conservation laws. A switch s is 
designed to turn V on (s = ON) or off (s = OFF). 

3. After a sufficiently long time Ti = ti — tq the switch 
is flipped (local quench Hog —)■ Hon)- 

4. At a later time T 2 the switch is flipped again (local 
quench Hon —> Hog), and so on. 

We point out that similar dynamics were already consid¬ 
ered e.g. in Ref.®®. 
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FIG. 1. Cartoon of the quench protocol. Full (empty) circles 
represent sites that are occupied (empty) at the initial time 
ro = 0 and the dot-dashed lines indicate an undefined number 
of sites. The switch s (light yellow triangle) decouples a site 
from the others. The (thick) red and (thin) green dashed lines 
represent connections when the switch is in the status of the 
corresponding color. The ellipse encircles the subsystem, A, 
where macroscopic observables, like (3), are defined. 



FIG. 2. Space-time picture of the light-cones propagating 
from the switch. Tags “ON” and “OFF” are placed at the 
times Ti when the switch is turned on and off, respectively 
(before ti the state is two-site shift invariant). The thick 
segment represents subsystem A at time t. The subsystem 
is split in parts belonging to different phases: in the green 
region EVs are approximately one-site shift invariant while 
in the orange region there is antiferromagnetic order. The 
crossover (orange to green) is still developing and the EVs of 
macroscopic observables vary almost linearly in time. 


a. Example. For the sake of clarity, we analyze one 
of the simplest dynamics that can be controlled by a lo¬ 
cal switch through the mechanism described above and 
explained more extensively in Sec. I. The reader can find 
further examples in Sec. IV. We consider a fermionic 
chain with Hamiltonian 

L 

Hxy = + ( 1 ) 

where satisfy the algebra {c.,Cj} = 5ij, {ci,Cj} = 0 
and periodic boundary conditions Cl+i = Ci are imposed. 


In (1) we omitted a multiplicative constant with the di¬ 
mensions of an energy, which can be viewed as implicitly 
attached to the time. We note that Eq. (1) is mapped to 
the XY spin-1 chain^^ by a Jordan-Wigner transforma¬ 
tion. It was shown^° that this model has local conserva¬ 
tion laws that break one-site shift invariance and are odd 
under a shift by one site. Clearly, they can exist only if 
L is even. Let us choose L even and set the Hamilto¬ 
nian iLoff = Hxy- The chain’s length can be effectively 
modified through a local perturbation. For example the 
following local term 

c|(C2 + CL+ 7(4 - cj)) - 4^2 - 7ci4 I ^ (2) 

decouples the first site and has the practical effect of 
reducing the chain’s length by one site. This in turn 
results in the extinction of infinitely many local conser¬ 
vation laws. We therefore set Hon = Hxy + V. Since V 
spoils only conservation laws that are odd under a shift 
by one site, we prepare the system in a state that breaks 
one-site shift invariance. In practice, this can be achieved 
by quenching from an antiferromagnetic or a dimerized 
phase. In the rest of the paper the pre-quench state will 
be the Neel state (occupied sites alternating with empty 
sites). We point out that in the thermodynamic limit 
the distinction between odd and even chains is reduced 
to the presence or absence of a defect in the initial state 
at the position of the switch. By analogy with the phe¬ 
nomenon of metastability in phase transitions, the defect 
acts like a seed crystal from which a larger crystal will 
grow. Fig. 1 sketches the system configuration. 

Globally different states can be distinguished by the 
EVs of macroscopic observables, i.e. local observables 
that are averaged over a sufficiently large number of sites. 
We focus on the staggered magnetization per unit length 
in a connected subsystem A 

= ^ (3) 

1^1 itA 

where |A| is the subsystem’s length, which we assume to 
be even. For even distances, at the initial time m^’^[A] = 
i. The time evolution of the Neel state under HqS can 
be easily worked out 


lim lim [HI 

|A|-).ooi-!-oo ^ 


1 ItI 
21 + 171 ■ 


( 4 ) 


On the other hand, time evolution under Hon results in a 
one-site shift invariant state (there are only one-site shift 
invariant local conservation laws), and hence 


lim lim = 0. (5) 

|A|—>-oo t—xx) 


Our protocol of non-equilibrium dynamics is more com¬ 
plicated than a global quench but these limits are still 
appropriate reference values. 
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FIG. 3. The staggered magnetization per unit length in the 
subsystem consisting of the sites i € [2,7] (|A| = 6, d = 0) 
for 7 = 2 and various chain’s lengths. The dotted vertical 
lines correspond to the times when the switch is activated 
(thin green) and deactivated (thick red); a faint green shadow 
highlights the intervals when the switch is on. The dotted 
horizontal lines are the reference values (4) and (5). The 
dashed black curve is the prediction (6). 

b. Analysis. The dynamics of macroscopic observ¬ 
ables in A are determined by four parameters: 

(a) The chain’s length L; 

(b) The subsystem’s length |A|; 

(c) The distance d between A and the switch; 

(d) The time Ti = n — n-i that elapses between consec¬ 
utive changes of the switch status. 

The role of the parameters can be understood using a 
semiclassical picture^^ based on the motion of quasi¬ 
particles excitations originated from/scattering with the 
switch. The quasiparticles travel at finite speed, bounded 
from above by the Lieb-Robinson velocity vm = HtI ~ 1|- 
Therefore the light-cone emerging from the switch prop¬ 
agates at Vm- Inside the light-cone a stationary state 
compatible with the perturbation is developing while, 
outside, the state time evolves as in the absence of the 
perturbation (c/. Fig. 2). Around the light-cone there is 
an unstable transient region. 

Our analysis shows that the finite size effects, (a), be¬ 
come visible after a time t ~ L/ (2vm) but the qualitative 
behavior is almost unchanged until a time t ^ L/vm (see 
Fig. 3). 

The main effect related to the subsystem’s length, (b), 
is in the time that the light-cone employes to flow through 
the subsystem: the larger A is and the slower the transi¬ 
tion between one state and the other (c/. Fig. 4). During 
this time the EVs of macroscopic observables experience 
a typical linear time evolution: the light-cone, moving at 
the Lieb-Robinson velocity, splits the subsystem in parts 
that approximately belong to different phases. 

The distance d from the switch, (c), introduces a de¬ 
lay St = djvM between the time at which the switch 
is flipped and the time at which observables in A start 



FIG. 4. The staggered magnetization per unit length for var¬ 
ious subsystem’s lengths starting from the second site {d = 0) 
with the same quench parameters and notations of Fig. 3 in 
a chain of 400 sites. The dashed lines are the predictions (6). 



FIG. 5. The staggered magnetization per unit length in the 
subsystem [2 -f d, 13 -I- d] with the same quench parameters 
and notations of Fig. 3 in a chain of 400 sites for various 
distances d. The dashed lines are the predictions (6). The 
time is shifted to compensate for the finite velocity at which 
information propagates. 

feeling the change. In addition, only the quasiparticles 
with velocity v > um/( 1 + have reached the 

subsystem, so the transition is also slower at larger dis¬ 
tances, where the effects of the unstable regions around 
the light-cones are amplified (see Fig. 5). 

The time Ti between consecutive local quenches, (d), 
puts a limit on the time left for the subsystem to respond 
to the change (c/. Fig. 6). 

Following a semi-quantitative reasoning based on the 
semiclassical picture of Fig. 2, we conjecture an asymp¬ 
totic expression for (see Sec. II), which is expected 
to hold in the thermodynamic limit when all the param¬ 
eters are sufficiently large: 

^ 1 l7l r dfc 7^sin^(fc/2)/x^^^^(fc) 

2 1-b ItI 27r 1 -b COS fc-b 7^(1 — cos fc) 

Here {r} is the set of the times Tj at which the switch 
is flipped; the details of the dynamics are encoded in a 
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FIG. 6. The staggered magnetization per unit length in the 
subsystem [2, 7] (|A| = 6, d = 0) with the same quench pa¬ 
rameters and notations of Fig. 3 in a chain of 400 sites for 
various time intervals Ti = T between successive switches. 
The dashed lines are the predictions (6). 

single function 

[min(2|4|(t-Tj),|A| -bd) - niin(2|4|(t - Tj), d)] (7) 

where £k = (cos^ | + 7^ sin^ is the dispersion rela¬ 

tion of the two species of quasiparticles that diagonalize 
HoS (1) with the momentum k defined in terms of trans¬ 
lations by multiples of two sites. The numerical data 
are in good agreement with (6), which indeed captures 
the relevant features of the dynamics in all the situations 
investigated (c/. Figs 3, 4, 5, 6). 

Although the quench details affect the time evolution 
of macroscopic observables (as a consequence of the tran¬ 
sitions being always incomplete), the entire process shows 
surprising reversible features: the profile of the staggered 
magnetization is only weakly influenced by the past ac¬ 
tions on the switch. This can be inferred also from (6) 
and (28): when t S> tj the contribution to from 

the switches at r, < rj is reduced to j which 

has the practical effect of resetting the EV to (4), for j 
even, or (5), for j odd. 

c. Projective measurement. Here we investigate 
whether our construction is stable under projective mea¬ 
surements (see Sec. III). Fig. 7 shows for |A| = 6 
following repeated measurements of all the occupations 
numbers Uj = cjcj in A. Generally, after the measure¬ 
ment there is a time interval characterized by sizable fluc¬ 
tuations that make the time evolution noisier. Neverthe¬ 
less the qualitative behavior of does not change. In 
average the measurements distinguish the two phases. 

d. Summary and discussion. We proposed a general 
protocol of non-equilibrium evolution that could be used 
to engineer a quantum device in which global properties 
of a state are controlled by a local switch. We tested 
the protocol on a quantum XY fermionic chain, where 
we analyzed the staggered magnetization averaged in a 


{ 0 , 1 . 0 , 1 , 1 . 1 } { 1 . 0 , 0 , 1 , 0 , 1 }{ 0 , 1 , 1 , 1 , 0 , 1 } { 0 , 1 . 0 , 1 , 1 . 1 } { 0 , 1 , 1 , 1 , 1 . 0 } { 0 , 0 , 0 , 1 , 0 , 1 } { 0 , 1 , 0 . 1 , 0 , 1 } 

{1.1.1,0,0.1} {1.1,0,1,0,0}{0,0,1,l,0.1}{0.i,1,1,1,1} {0,1,0,0,0.1} {0,0,0,1.1,1}{0,1,1.0,1,1} 



FIG. 7. The staggered magnetization per unit length in the 
subsystem [2,7] (|A| = 6, d = 0) with the same quench pa¬ 
rameters and notations of Fig. 3 in a chain of 800 sites. The 
squares (red, when s = OFF, and green, otherwise) represent 
(simultaneous) projective measurements of all the occupation 
numbers in A, repeated at time intervals of 24 starting from 
t = 18 (the outcomes are shown on top with the color of the 
squares). The dashed curve Is (6) without the measurements. 


subsystem and conjectured an asymptotic expression for 
its time evolution. In order to assess the stability of the 
setup we kept all the parameters finite and tested the 
protocol under projective measurements. 

The effect described relies on the absence of a one-to- 
one correspondence between (the bulk part of) the local 
conservation laws with and without a local perturbation. 
These situations are very common and emerge also in 
generic models with a global symmetry {e.g. particle 
number conservation) if the switch is associated with a 
symmetry-breaking defect (see Sec. IV for an explicit ex¬ 
ample). The generality of this situation opens the door 
to experimental realizations. 

In integrable models the switch could change the 
boundary conditions from periodic to open, where the set 
of local charges is almost halved^^’^® (see also Sec. IV). 
More general local perturbations would spoil integrabil- 
ity, as could be inferred also from Ref.^®. 

Finally, we note some similarities with the protocol of 
joining together two chains to form a single one (and 
other similar dynamics), considered e.g. in Ref.^^. Also 
in that case a light-cone that separates globally different 
states is originated from a given point (the junction). 
The two situations contrast in the initial states, which 
we assumed to be homogeneous and non-stationary in 
the bulk. 
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SUPPLEMENTAL MATERIAL 

I. LOCALIZED DEFECTS IN SPIN CHAIN HAMILTONIANS 

Let us indicate with Hq a generic homogeneous Hamiltonian with periodic boundary conditions and with Hi the 
same Hamiltonian with a defect d(ro) localized around position rg: 

Hi=Ho+ d{ro). (8) 

We consider homogenous initial states Idio) and focus on situations where the expectation values of local observables 
relax to stationary values both under Hq and under Hi. 

An important consequence of local relaxation is that the expectation value of the commutator between the Hamil- 
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tonian and any operator O with a finite typical range approaches zero in the limit of inhnite time, indeed 

dt ^ 0 ^ ^0. (9) 

Let us now consider a conservation law Q of Hq that can be deformed close to the defect so as to be conserved also 
for s = 1 : 


[Q,iLo]=0 [Q + SQ,Hi]=0 ||dQ||<oo, (10) 

where || • || is the operator norm (the maximal eigenvalue in absolute value). 

The variation of the time evolution of Q under Hi is given by 

I (vI/o|e*^^‘Qe-*"i‘ - QlM-o) I = I (^ol^Q - I < 2 || dQ |i . (11) 

Being the norm of 6Q finite (see (10)), the expectation value of Q, which is generally proportional to L, remains close 
to the original value at any time after the quench. 

In the physical situations when the infinite time limit of the expectation values of local observables commutes with 
the limit of large distance from the defect, it is reasonable to expect that the stationary behavior of observables can be 
described by the statistical ensemble that maximizes the entropy under the constraints of the integrals of motion (this 
is indeed the principle that holds in the limit limi_,.oo limd_>oo). In all the cases studied we found that this happens 
when all the conservation laws that are relevant to the stationary properties for s = 0 are simply deformed by the 
defect 




SQj ■ 


( 12 ) 


On the other hand, the situation becomes much more involved if there are charges of Hq that become extinct. For 
example, let us assume that there is a charge Q of Hq such that 


lim (4'o|e*-^i‘i[d(ro),Q]e 7 ^ 0 

^ t —^00 


(13) 


Since 


[Hi,Q] = [Hq + d(ro), Q] = [d(ro), Q ], 
the variation of the expectation value of Q is given by 

I - Q|vI/o) | = | /‘ dr (vI/o|e*^^"*Kro), Q]e-*^^n«'o) I • 

The right hand side can be bounded from below as follows: 


(14) 


(15) 


dT(4'o|e*^^'^i[d(ro),Q]e 


dr (vI/o|e*^^"i[d(ro),Q]e-*^^n'f'o) + dr (vI/o|e*^^"z[d(ro),Q]e-*^^n'f'o) 


> -T II z[d(ro), g] II +(t - T) min | {M>\A^^A[dirQ), Q]e-^^^^\^) \, (16) 


where T is sufficiently large that 


min I (^o|e*"^"*Kr-o), Q]e-*"^"|^o) I > 0 . (17) 

t>T 

Since || i[(i(ro), Q] || is finite (the operator is quasilocalized around the defect), this gives a lower bound to the variation 
of the expectation value of Q that increases with the time, diverging in the limit of infinite time. By (11), there can 
not be a corrective term SQ that satisfies (10), and hence there is no quasilocal deformation of Q that makes the 
operator commuting with Hi: Q is destroyed by the defect. A linear increase (or decrease) of the expectation value 
of a charge of Hq turns out to be the typical signature of its extinction. 
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We now show that a nonzero Aq, (13), results in the propagation of a light-cone that separates regions with globally 

different properties. To that aim we consider the density of the extinct charge Q, where index I signifies that the 
operator acts nontrivially only around site I 

Q = ^qi- ( 18 ) 

£ 

Since dt (^'d^l'I't) = summing over the sites in an interval [d-|-T q, n-|-T q], with 0 < d < n, gives the 

continuity equation 

n+ro n+ro 

^ QlW'i’t) = {'i’t\Jn+ro - Jd+rol'i't) ■ ( 19 ) 

£—d-\-rQ i—d-\-rQ 

Here we used that the extinct charge is conserved in the bulk, and hence the commutator can only depend on surface 
terms, which we called J (they are local operators). We note that, as long as d is larger than the sum of the ranges 
of Hq and Q, the definition of J does not change if Hi is replaced by Hq. Let us now consider the limit d vt n, 
with V the Lieb-Robinson velocity. We find 

d<\i—ro\<n \i—rQ\<d |£—ro|>n 

d^{^^,,\^[d{ro),Q]\'ft)-^^{'^'t\ Y mt)^AQ. (20) 

|^-T-o|<d 

In the first step of the second line we used that the density is practically constant until the information about the 
defect has not yet reached the site; in the second step we instead used that, at times t ^ d/v, local observables relax. 
Putting all together we get 

(«'t|(J„+,„ - Jd+ro) + {Jro-d - dro-n)|«'t) Aq 0 , (21) 

where, we stress, Jn+ro and Jro-n {Jd+ro and Jro-d) are outside (inside) the light-cone. This means that the 
expectation value of Ji depends on whether ^ is inside or outside the light-cone, even for d^ 1. 

One can use similar arguments to show the effects of flipping the switch. For example, let us assume to turn off the 
switch at time T 2 . Because of the Lieb-Robinson bounds, at times t < +2 -\- d/v the expectation value of Jd+ro (and, 
in turn, of Jn+ro) ‘^an not be changed. On the other hand, in the opposite limit d ^ n v{t — T 2 ) + d ^ vt, it is 
reasonable to expect relaxation to the values corresponding to time evolution under Hq, which in fact are the initial 
values. Let us then consider the remaining two nontrial limits 

1. d < v{t — T 2 ) + d vt n; 

2. d <C v{t — T 2 ) + d n vt. 

In case 1 we find 

d<\£—ro\<n \£.—rQ\<d 

where we used that now Q is conserved and hence (5't|i[d(ro), Q]|'I't) —>■ 0. This results in 

{'i't\{Ju+ro - Jd+ro) + {Jro-d ” ^-n)|^t) 0 . (23) 

Since {Jn+ro) ^ad still the initial value, this suggests that also Jd+ro returned to the original value. 

Case 2 can be deduced from the other limits: Jn+ro turns out to have the value corresponding to Hi, whereas Jd+ro 
already relaxed to the initial value. 

We stress that we did not use any specific property of the Hamiltonian, so this physical picture applies to nonin¬ 
teracting, interacting integrable and even generic models with at least a local charge besides the Hamiltonian. 

In conclusion, in the presence of extinct charges we should expect that the stationary state that describes the late 
time properties under Hi deep inside the light-cone is macroscopically different from the stationary state that gives 
the correct behavior outside the light-cone (which can be simply obtained replacing Hi by Hq). In other words, the 
limit of infinite time does not commute with the limit of infinite distance. In addition, flipping the switch corresponds 
to reverse the situation. 








II. SEMICLASSICAL FORMULA FOR (SIMPLE) MACROSCOPIC OBSERVABLES 


Let us focus on the model (1). The space-time picture of Fig. 2 can be easily translated into a semiclassical formula 
for the time evolution of macroscopic observables. We consider extensive quadratic operators 0~ that are odd under 
a shift by one site and restrict them to a subsystem A as follows 


^ 1 tr^[e> ] 0 I^ 
1^1 tr^pT] 


(24) 


where A is the complement of A. The staggered magnetization per unit length considered in the main text 
belongs to this class. 

If the time ri of the first switch is sufficiently large, we can assume {O'^) to have approximately reached the 
GGE value associated with time evolution under Hog; this is an integral over the momentum (defined in terms of 
translations by multiples of two sites) of a function depending on the observable 



lim (4-01 

>-oo 


O^e- 



dk 

27T 


f{k;0-). 


(25) 


We interpret f{k; O ) as the contribution given by the quasiparticle with momentum k to the expectation value. In 
particular, for the staggered magnetization we find 




7 ^ sin^(A:/2) 

1 -I- cos/c -I- 7^(1 — cosfc) 


(26) 


On the other hand, under Hon the expectation values of odd observables approach zero. 

The first time that the switch is flipped (s —)■ ON) a destructive light-cone emerges. This is responsible for a 
negative contribution —f{k;0~) when the quasiparticle is (semiclassically) inside A. Flipping again the switch the 
new light-cone is instead associated with positive contributions f{k;0~), and so on. This results in 





fj{k;0-) 


1 + kf,{T}{k) 


(27) 


where given by 

k-t{r}{k) = - Tj)[min{2\e'^\{t - Tj),\A\ + d) - min(2|4|(t - Tj), d)] , (28) 

counts the number of quasiparticles with momentum k in A (per unit length) with an overall sign depending on the 
status of the switch when the quasiparticles were produced. 

For quadratic operators that are even under a shift by one site, the semiclassical approach does not show any 
effect caused by the switch. We do not expect that this holds true for non-quadratic (even) operators, since the 
decomposition of their expectation values (EVs) by Wick’s theorem involves also EVs of operators that are odd under 
a shift by one site. 


III. PROJECTIVE MEASUREMENTS 


A (von Neumann) projective measurement n[0] of an observable O involves the collapse of the wave function 


0 = J2KPx, {Px^Px, 

i 




^ijPXi: f — j) 

Pi^i[0] = K) = i^lPx-J'f) , 


(29) 


where p{p[0] = A) is the probability of measuring A. In principle, the collapse could disrupt the nice local prop¬ 
erties that characterize our protocol of non-equilibrium time evolution. The response of the system to projective 
measurements is therefore a severe test for the dynamics considered in this paper. 

Since the Hamiltonian is quadratic and the initial state is a Slater determinant, we can restrict ourselves to pre¬ 
measurement states |4') that are Slater determinants. However, after a measurement the state becomes generally 
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very complicated anyway. For the sake of simplicity we consider measurements that preserve gaussianity, a sufficient 
condition being the applicability of Wick’s theorem to the density matrix for any j. A simple observable 

with this property is the occupation number Uj = CjCj. Its measurement results in the following collapse of the density 
matrix 


1^) (^1^ 


nj 1^) 

(l-nj-)l^) (^1(1-nj) 
1 - 


(p= mK] = 1 

(p= 1- fi[nj]=0. 


(30) 


Any expectation value can be written in terms of the two-point correlation functions, which are conveniently written 
in terms of the Majorana fermions ({ 0 ^, 0 ^} = 2Sij) 


02i-l — c\ + o, 


021 


= *{c* - 4) ■ 


(31) 


We call r the correlation matrix 


Fij — Sij ('kjaiOjl^') . 


(32) 


The probability of measuring ulrij] = with s = ±1, is p = ^ After the measurement the correlation 

matrix transforms as follows 


2/i[nj] = l+s 

t in -' 


1 - 


(<5 






1 - isT 


2i,2j-l 


+ 


. S2j/52j-l,n — <^2j-l.f<^2j,n + r„_2j-ir^,2i — r^^2j-ir„,2j ^oo^ 

----— 1 —--- • 

1 — ^si 2j,2j-l 


Since = 0, all the occupation numbers in a subsystem A {j € A) can be measured simultaneously (this is 

equivalent to measuring an observable with incommensurate aj). The resulting correlation matrix can 

be obtained applying several times the transformation (33). The probability of the outcome is obtained analogously. 
Fig. 7 shows the time evolution of the staggered magnetization for a particular set of results of the measurements, 
generated randomly with the respective probabilities. The reversibility properties of the dynamics along with the 
regularity at which the measurements are performed and the switch is flipped can be used to reconstruct a probability 
distribution for the staggered magnetization. As far as expectation values are concerned, our preliminary analysis 
suggests that it is as if the quantum state were almost perfectly replicated before each measurement. In particular, 
the two phases can be distinguished by the time average of the outcomes. 

Finally, we note that a measurement of n-r+i — Ur breaks gaussianity when the result is 0, indeed the projector on 
the corresponding eigenspace, i.e. 1 — rir — n^+i + 2nrnj.+i, is not gaussian (in terms of the fermions (31)). 


IV. FURTHER EXAMPLES 

In this section we show some numerical results for the following examples: 

1. Switch of the boundary conditions in a fermionic XY chain in a transverse field; 

2. Switch associated with a localized magnetic field that breaks a f7(l) symmetry in a nonintegrable spin-^ model. 


A. Case 1 


Let us include a magnetic field h to the Hamiltonian (1) 

A.. . A. 

e+u ■ 2 




X1[^(c^4-hi-4c^-u) + |(c^c^-U-44+i) + ^(24c^-1) -^\{clc\-c\ci) + '^{clCi-c\c\) , (34) 




where we imposed periodic boundary conditions Cl+i = ci. The field destroys all the extra (non-abelian) charges that 
characterize the XY model in zero field when L is even. The switch changes the boundary conditions from periodic 
(s = 0) to open (s = 1). 

As initial state we take the fermion vacuum 


qI^o) =0. 


(35) 
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FIG. 8 . The energy density at position n and time t corresponding to the time evolution of the state (35) under the XY 
Hamiltonian (34) with h = 0.333 and 7 = 1.732 (left) or 7 = 1.1 (right) in a chain of 300 spins (notice that the scale 
of the colors is the same). Left. There are extinct local reflection symmetric charges. In the limit 1 vt n the energy 
density is undistinguishable from the one corresponding to time evolution with Hq, which is also the initial value (black); in 
the limit 1 <C n <C ft a new stationary behavior emerges. At times t > 250 finite size effects become visible. Right. All the 
local reflection symmetric charges survive the boundary. The effect of the boundary fades away at large time and distance: the 
switch can not be used to control global properties. 



FIG. 9. The energy density at position n and time t corresponding to the time evolution of the state (35) under the XY 
Hamiltonian (34) with h — 0.333 and 7 = 1.732 in a chain of 300 spins. The switch is ON (s = 1) until t = 100, then it 
is turned OFF (s = 0). In the right panel the switch is turned on again at t = 200. Left. After flipping the switch the original 
value of the energy density is slowly restored. Right. After flipping the switch the second time a new light-cone emerges with 
the same morphology of the original one. 


We set h = 0.333 and consider two values for the anisotropy: 7 = 1.732 and 7 = 1.1. 

Generally, in integrable models, changing the boundary conditions from periodic to open destroys the charges that 
are odd under chain inversion^^. However, the initial state is reflection symmetric, so the extinction of the odd 
charges should not be relevant to time evolution. This example is still significant because Ref.^® showed that the 
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FIG. 10. The cumulative distribution function F{s) of the nearest neighbor spacing distribution {i.e. the probability density of 
the difference between two consecutive levels of the unfolded spectrum of the energy) for the Hamiltonian Hq (36) with A = 0.5, 
g = —0.5, U = 0.25 in the reflection symmetric sector with zero momentum, zero magnetization and spin flip erf —>■ 1 for 

■KS^ 

L — 20. The data are very well described by a Wigner distribution {F{s) — 1 — e ^), which is typical of nonintegrable models 
(in integrable models the typical distribution is instead Poisson F{s) = 1 — e”®). 

set of (quasi)local conservation laws for s = 1 depends on the ratio q — \h/{^^ ~ 1)1; specifically, for g > 1 all the 
reflection symmetric charges of the periodic chain survive the boundary, for q < 1 instead there are less conservation 
laws. The two values of 7 that we consider are representative of the two situations. Consequently, for 7 = 1.732 the 
switch is expected to control global properties, while for 7 = 1.1 it is not. 

Fig. 8 shows the main difference between the case q < 1 and g > 1. In the former, time evolution under Hq results in 
a stationary state different from the one with periodic boundary conditions (s = 0). On the other hand, for q > 1 the 
limit of large time seems to commute with the limit of large distance. This is perfectly consistent with the theoretical 
picture presented in Sec. I. 

In Fig. 9 we check reversibility by turning off the switch and flipping it again after a sufficiently large time. The 
novel light-cone has the same morphology of the original one. 

We would like to point out that the particular profile of the energy density is probably model dependent and, for 
example, the rather complicated behavior in the limit of infinite time at fixed ratio j could be a consequence of 
integrability or of the absence of interactions. Nevertheless, as shown in Sec. I, the noncommutation of the large time 
limit with the large distance turns out to be a general property of any model with a defect that destroys a relevant 
set of conservation laws. 


B. Case 2 


As a final example, we consider a nonintegrable model with Hamiltonian 


Hs 




e=i 


+ -I- 5cr|cr|^2 + U{(T( 


) —StT 


X 

1 5 


(36) 


where periodic boundary conditions = uf are imposed and s takes the value 1 or 0 depending on whether the 
switch is ON or OFF. We set the Hamiltonian parameters as follows: A = 0.5, g = —0.5 and U = 0.25. Fig. 10 shows 
that this choice makes the model generic. For s = 0 the Hamiltonian is invariant under rotations about z, indeed it 
commutes with | ^ 1 - 

In the following we provide numerical evidence that the local magnetic field in the x direction, which breaks 
rotational symmetry around z, is responsible for the emergence of a front that separates two regions with different 
magnetization. 
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FIG. 11. The variation in the total magnetization | (ct|) after a quench from the state (37) under Hamiltonian (36) 

with A = 0.5, g = —0.5, U = 0.25 and s = 1, for several chain sizes. The vertical dashed line is an estimate of the time at 
which finite size effects become significant even for L — 14. Data are consistent with a linear decrease of M^. In the inset, the 
local magnetization = | {o'f}, which is proportional to the time derivative of M^. 
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FIG. 12. The profile of the local magnetization | (a^) after a quench from the state (37) under Hamiltonian (36) with 

A = 0.5, g = —0.5, U = 0.25 for the chain sizes L = 8,10,12,14. In the upper panels the switch is always ON (s = 1); in the 
lower panels the switch is turned off at t = 1.3. 


As initial state we take the product state 


l^-o) = It •• • t) (37) 

We study time evolution by exact diagonalization techniques, which allow us to access only small systems (we inves¬ 
tigated chains up to L = 14). 
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FIG. 13. The profile of the local magnetization = 1 (cr^) after a quench from the state (37) under Hamiltonian (38) with 
A = 0.5, g = —0.5, U = 0.25 for the chain sizes L = 8,10,12,14. In the upper panels the switch is always ON (s = 1); in the 
lower panels the switch is turned off at t = 1.3. 


The total magnetization does not commute with Hi. A generic defect that does not commute with is expected 
to destroy the conservation law . As discussed in Sec. I, the standard signature would be the presence of a nonzero 
slope in the time evolution of the total magnetization in the z direction. Although the hnite size effects are significant, 
Fig. 11 is compatible with a linear decrease in time (which should be regarded as the standard situation) and hence 
we deduce that there is no term (quasi)localized around the defect that makes commuting with Hi. 

The upper panels of Fig. 12 show the profile of the local magnetization m| = ^ {af) as a function of the time for 
various system sizes (to give some idea of the finite size effects). One can fairly notice the emergence of a macroscopic 
region with a different local magnetization. The lower panels correspond to turning off the switch at the time 
t = 1.3. One can perceive the propagation of a new light-cone that is going to restore the original magnetization. The 
dimensions of the systems investigated are clearly too small to deduce any behavior purely from numerics, however 
the data are consistent with the development of regions with globally different properties, as explained in Sec. I. 

Fig. 13 shows the results of the same analysis for a different switch that turns on a local magnetic field in the z 
direction 


H2= Ho-at. (38) 

This time the Hamiltonian remains invariant under rotations about z and we do not expect macroscopic effects. 
Fig. 13 is compatible with our belief, in the sense that the behavior inside the light-cone appears closer to the one for 
s = 0 (the initial value). 














